Zeta function regularization 
in de Sitter space: 
the Minkowski limit 



Alan Chodos^ and Andras Kaiser^ 
Center for Theoretical Physics, Sloane Physics Laboratory 
Yale University, New Haven, Connecticut 06520-8120 



^ chodos@yalph2.physics.yale.edu 
^ kaiser@hep Unix. physics . yale . edu 



Abstract: 

Wc study an integral representation for the zeta function of the one-loop effective potential for a min- 
imally coupled massive scalar field in D-dimensional de Sitter spacetime. By deforming the contour of 
integration we present it in a form suitable for letting the de Sitter radius tend to infinity, and we demon- 
strate explicitly for the case D=2 that the quantities C(0) and C'(0) have the appropriate Minkowski limits. 



I. Introduction 

In this paper we examine some properties of the effective potential for a minimally coupled scalar field 
in de Sitter spacetime. We regard this as a preliminary study with application to symmetry breaking effects 
in de Sitter space, which can have significant cosmological consequences [1]. These can be investigated if one 
also includes a potential for the scalar field. 

Our particular concern in this work is to examine the rcgularization of the effective potential by the 
zeta- function method [2] , and to show that after such rcgularization one can compute explicitly the limit of 
the potential as the de Sitter radius a tends to infinity, and recover thereby the Minkowski space result. 

Of course de Sitter space itself becomes Minkowski space as the radius becomes infinite, and furthermore 
at the naive level the limit of the de Sitter zeta function C(s) goes over into the Minkowski one (see below). 
However in order to regularize the effective potential one must perform an analytic continuation of C to a 
neighborhood of s = 0. After this this is done, it is no longer evident (or even necessarily true) that the 
zeta- function possesses the correct a ^ oo limit. 

In other words, the process of analytic continuation and the a ^ oo limit may not commute with each 
other, and it is reassuring to check explicitly that the analytically-continued C(s) still possesses the expected 
Minkowski limit. 

A simple example of what needs to be checked is provided by the series 

oo 

/(s; a) = ^ e-("+i)"(*-2)(-l)" (1.1) 

n=0 

(a > 0), which for 3?s > 2 clearly obeys 



lim /(s;a) =0 (1.2) 

a— >oo 

But /(s; a) possesses an analytic continuation in s 

f(s;a) = (1-3) 

which, in a neighborhood of s = has the limit 



lim/(s;a) = l (1.4) 

a—*oo 

Let's take a look at the Minkowski zeta-function first. In D dimensions, it is defined by 

*) = Ejer^ (1-5) 

where /i is a scale factor with dimensions of mass, that is included to make the dimension of ( independent 

of s. 

In Minkowski spacetime the sum is really continuous, so we can replace J2k ^^^^ VdPk/{2'K)^ {V is 
the volume of the spacetime) to transform it into an integral form 

or divide by V to obtain the C density 
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^^'^ ~ J (27r)^ • ^^■^> 

II. de Sitter zeta-function 

Proceeding to the de Sitter case we note that D- dimensional de Sitter spacetime can be cinibciddcxi in 
a D+1 dimensional Minkowski-space (we use signature < +,—,—,...,— >), where it is a D dimensional 
hypersurface determined by 



^o2 



E^^' = (2.1) 



We use the following parametrization of this space 



zq = sinh(a;o) 

zi = cosh(a;o) cos(a;i) 

Z2 = cosh(a;o) sin(a;i) cos(a;2) 



(2.2) 



We start with the de Sitter action 

S = J d^x !^ 5^0*9"^ - mVl^ I (2-3) 
Substituting xq = —ix4, we obtain the Euclidean action 

S = J d^x I d^4>*d^(t> + m^\<f>\'' I (2.4) 
The classical (j){x) field configurations form a vector space with the scalar product 

<(f>,il)> = j dx (l)*{x)ip{x) (2.5) 



and the second variation of the action 

5^S 



6(l){x)6(j)* {y) 



= {De W) Six - y) (2.6) 



is similarly considered a matrix (bilinear operator on the vector space). The zeta function is defined 
with the eigenvalues of this matrix 

C(^) = (2-7) 

i ' 

where /x^* is included to get a dimensionless quantity. When we substitute xq = —ixi, cosh(a;o) trans- 
forms into cos(a;4), and sinh(a;o) into — isin(a;4). This means that now we have 

Zq = —i sin(x4) 

Zi = cos(a;4) cos(a;i) 

(2.8) 

Z2 = cos(a;4) sin(a:i) cos(a:2) 



which shows that Euclidean de Sitter space is a (series of) D-sphere(s), therefore the eigenfunctions 
(eigenvectors) of 5'^S/5(j){x)5(j){y) are generalized spherical functions, with the eigenvalue 
\ = l{l + D — V)lo? + rr?. The multiplicity of the modes with the same / is [3] 
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r{D + i + i) 



T{D + l-l) 



T{D + 1)T{1 + 1) r{D + l)T{l-l) 



so 



as) = E 



T{D + l + l) 



T{D + l-l) 



T{D + 1)T{1 + 1) T{D + 1)T{1-1) 



2s 



l{l+D-l) 



(2.9) 



(2.10) 



T{D + 1) 



where 



^ E [(D + ;)(D + Z-1) - Z(l-l) 



r{D + i-i) 1 

r(/ + i) [ (/ + ^)2 - /32 



(2.11) 



Let's make a brief digression here, and check how the ^-density approaches the Minkowski case as a — » oo. 
Divide ( by the volume of de Sitter space, a^Vo, where Vd is the surface area of a i3 + 1-dimensional sphere, 
in other words the volume of the manifold . 



a^Vo T{D + 1) 



^ J2 [{D + l)iD + l-l)~ 1(1-1) 



TjD + l-l) 
T{1 + 1) 



,2s 



(2.12) 



Transform the sum over discrete values of I into a continuous integral by substituting l/a^ k and 1/a 
dk. We also assume / > 1, yielding (£)+/)(£)+/- 1) - l{l-l) 2DI = 2Dka,T{D+l-l)/T{l + l) Z^-^, 
and l{l + D - l)/a'^ ^ P. This gives 



1 



(is) 



2D 



T{D + 1)VdJ (k 



„D-l 



dk 



2D ^2 



T{D + 1) Vd Vd-1 J (P + m2) 



(2.13) 



Now substitute the value Vd = 27r— /T{^) 



1 ^ r(f ) r(^) 



^2s p(D±2) r(i?±i) r d^k 



u^Vd'^'"' r{D + l) in^+y^ J (fc2+TO2)s T{D + 1)ttD+^/^ J (fc2 + m2) 

And finally use the identity r(^)r(^ + l/2)/r(2^) = 2i-2^7rV2 

rf^fc r d^k /x2« 



2s o-D 



1 /;2« 2 



(F +7^2)^* 



{27t)d (A:2+m2)« 



which is the Minkowski (^-density. 

Now back to the de Sitter zeta-function, use the transformation [4] 



[a^-P^r r(7)7 V2/3 



e-"* /^-i/2 (/?i) rft 



(2.14) 



(2.15) 



(2.16) 



We get 



oc —1/2 
' 

Concentrate on the terms involving I, and evaluate the sum 
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^ T{D + I - 1) ^_i(o+2;-i)t 



= -2D 



dt 



= -2D 



dt 



T{D-l) 



t\D-l 



= T{D + 1) 



2cosh(t/2) 



[2sinh(t/2)]^ 

Now use this result in the expression of the zeta- function 



(2.18) 
(2.19) 
(2.20) 
(2.21) 



C{s) = 



T{D + l) T{s) 



t 

2^ 



^ [2sinh(V2)]^ 



{(3t) dt 



(2.22) 



With an appropriate correction factor, one can extend the contour of integration from —oo to oo, 
integrating above the origin. We will need to introduce a branch cut somewhere in the complex plane 
between the origin and infinity; let this branch cut be on the negative imaginary axis (i.e. from the origin to 
—ioo), which corresponds to the contour passing the origin on the positive imaginary side to avoid the cut. 



C(5) 



1 



oo+'iA 



r(.) 



\ — ^in{2s — D) 



-cso+iA 



t 



s-l/2 



2cosh(t/2) 
[2sinh(t/2)]^ 



/2 m dt 



(2.23) 



Here A is an arbitrary positive quantity, but less than 27r to stay between the origin and the first pole 
of the integrand on the positive imaginary axis. This expression is the desired analytic continuation that is 
well-defined in a neighborhood of s = 0. 



III. Deforming the contour 

Use an integral representation of the Bessel-function [5] 



C(«) 



1 



00 + 2 A 



r(s) 1 - e''^(2^-D) 



-oo+iA 



t 

2p 



s-l/2 



2cosh(t/2) 
[2sinh(t/2)]^ 



(3t 



-1/2 I g/3t« 
J (1-«2)1 



-du 



dt 



(3.1) 



cso+zA 



[r(s)]2 l_e"(2«-i5) J \2 

— oo+zA 



^ 2cosh(t/2) 



[2sinh(t/2)] 



D 



O0tu 



(1-U2)1 



du dt 



(3.2) 



Wc want to expand this expression in s around s = 0, so that wo can deduce C and its first derivative 
from the result. These arc the quantities that determine the eflFective potential. 
First examine the integral over u: 



(1 - U2) 



2^1-s 



du — 



(l-u2)l-s 



du 



p/3t(l-«) _^ ^-0t{l-v) 

1)1-^(2 - «)!-« 



dv 



(3.3) 
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But 



J + if^t -0t ) (3.4) 



1 . 1 1 

e 

~ ' ~ ~ ' 

The first term yields 

1 



The second term can be written as 

1 



+ ln2 + 0(s) (3.6) 



1 ^ {-I3tv) 



^ ' n=l 



[l + .ln. + 0(.^)] [l + «ln(2-.)+0(.^)] 



dv (3.^ 



Since we are looking only for the first two powers of s in the expansion, and we have a 1 / 2s term 
already, we will not need the terms proportional to s. Therefore we only need 

/^E^^- (3-9) 

"-^ 

which gives 

E^O-'-^-i:^ ,3.0) 

n=l ^ m=l ^ 



So 



oc+'iA ^ - 

/x^^Q^" 1 /• /A 2cosh(^/2) 

~ [r(s)]2 l-eM2a-i?) y 1^2^ [2sinh(i/2)]^ 



-oo+iA 



p/3t 



(-/3i)' 



n=l 



ln2 



m— 1 



+ 0(a) 



+ (/3i 



^ (3.11) 



Here we have not yet expanded the terms outside the curly brackets in powers of s, we shall do it later. 
Evaluate the above expression for ma » 1. This means (3 = ^J{D — 1)^/4 — rri^a^ k. ±ima. 

Let's find a way of closing the contour of integration at infinity in the upper or lower half plane, as 
appropriate for the different terms. First analyze the expression 



Amat 



1 

2s 



In 2 



E 

n=l 



(—imat)"' 



-)n— 1 



In 2 



E 

m=l 



■yn—m—1 



m 



0{s) 



(3.12) 



We can write In 2 as X^^=i ^ "^/m, and with this. 
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E 

n=l 



{—imaty 



n-l 



~\n—l 



m=l 



m 



E 



{—imaty 



E 



m 



(3.13) 



E°° {—imat)"' ^ 2 
, (n + 1)! 1 

n=l ^ ' 



m=l 



+ 



m-2 
n+1 



(3.14) 



„imat 



which clearly shows an e / {—imat) asymptotic behavior for \mat\ 1. This combined with the 



in front yields l/{—im,at), which cancels with the similar term from jSt 



'l3t, so we can subtract 



it from the expression. Then the remaining expression will go to zero at least as fast as l/{m^a^t^). This 
allows as to close the contour for this term in the upper half plane. The integrand has poles in the upper 
half plane at 27rm, n= 1, 2... . Applying the residue theorem, we obtain the integral as a sum 



oo+«A 



[r(s)]2 1 - e»'^(2«-i5) 



/ 



-oo+zA 



2s-l 



2cosh(t/2) 
[2sinh(i/2)]^ 



-imat 



—imat 



+ — + ln2 + y 
2s 



{—imaty 



^n— 1 



In 2 



n-l 

E 

m=l 



2"- 



m 



(3.15) 



+ 0{s) 



dt 



[r(s)]2 i-ei^i^^-D) ^ r{D) dt^-^ 



E 



1 d 



D-l 



t=27rin 




2s-l 



2cosh(t/2) 



t — 2iTin 
2sinh(t/2), 



D 



(3.16) 



AmoLt 



-imat 



1 . . {-imatY ( , . ^ 2"-™-i 



—imat 2s 



n=l 



2"-Mn2 - E 



m 



0(s) 



Prom this form it can be seen that this contribution is finite, and it goes to zero with ma oo. 

For the e"*™"* term, we could similarly close the contour in the lower half plane, except that the 
integrand has a cut from the origin to —ioo, because of the {t/2y^^^. So the curve that we use to close the 
contour must come back from —ioo to in the positive real half plane, then go around the singularity at the 
origin anticlockwise, and go back to —ioo in the negative real half plane, before continuing the arc to — oo. 
The complete contour constructed this way does not encircle any singularities, therefore the loop integral is 
zero, which means, that the contribution of the piece of contour along the cut cancels the contribution of 
the piece above the real axis, which is the integral that wc want to determine (the two quarter-arcs from 
oo to —ioo and from —ioo to — oo contribute zero, similarly to the semi-arc in the upper half plane in the 
previous case). So our integral above the real axis will equal to the integral from —ioo to —ie on the branch 
continuous with the left side of the cut, then a circle around zero clockwise with radius e, and then back 
from —ie to —ioo. 

Now back to the previous expression, let's expand it in powers of s. 



J(s) 



oc+iA 



[r(s)]2 i-eM^s-D) J y2 

— oo+zA 



t\^" ^ 2cosh(V2) 



[2sinh(t/2)]^ 



(3.17) 



imat 2s 



+ + ln2 



E 

n=l 



{im,at)" 



-)n— 1 



m 



0(s) 



dt 
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-2iTi 



1 [s + s- 



'(21n(Ma) -2r'(l)+i7r) + 0(s3)] J 0[1 + 2sln^ + 0(6 



2cosh(t/2) 
[2sinh(V2)]^ 



(3.18) 



-zmat 



imoff 2s 



+ _L + in2 + y 2"-iln2 - V 

2s n! V 



m=l 



m 



+ 0(s) 



IV. Evaluation of the integral 

Let's split the integral into the following terms: 

/• 2 2cosh(V2) , 1_ 

' J t [2sinh(i/2)]^ 2s 

C 



2o.^^ 2cosh(V2) j_ 
t 2 [2sinh(i/2)]^ 2s 



2 2cosh(t/2) -e'""* 
i [2sinh(i/2)]^ imat 



U=l -T#^^?!^e---ln2c^t 
' t [2sinh(V2)]^ 



(4.1) 



/ - '^"rf/l e--- V f 2"-ln2 _ g ^ ) 

y i 2sinh(V2)P ^ n! V m / 

\ / /J n=l ^ m=l ^ 



keeping only the two lowest-order terms in s. Only I2 has a cut along the negative imaginary axis, so 
the other terms will decompose into the contributions of a series of poles along the negative imaginary axis. 
The contribution from all the poles except the origin can be expressed as 



[r(s)]2 1 _ ei<-^^-D) r(D) 



1 d 



D-l 



t= — 27rin 



2s-l 



2cosh(f/2) 



2sinh(t/2) ) 



D 



(4.2) 



-imat 



oimat 



imat 2s 



n=l 



n! 



2"-Mn2 - 



m 



0(s) 



which is finite and goes to zero with ma — > cxd, similarly to the contribution of the poles in the upper-half 
plane. 

Now calculate the conribution from the pole at the origin for all but 



2s j 



2 2cOsh(V2) ^-imat^^ 



2sJ t [2sinh(i/2)]^ 



hi' 



2s I ^ 



1 + 



1 D 
8 ~ 24 



1 — imat — -m^a^t^ + 



(4.3) 



- {-2m) ( \m^o? 



= 2m ( m^a^ - 1/12 ) 



/3 



/ t 



2 2cosh(t/2) -1 



t [2sinh(i/2)]^ imat 



dt 



/ dt -r 

ima I 



^ + ^'--24 



(4.4) 



Ia = ln2 



/ 



2 2cosh(V2) 
t [2sinh(i/2)]'° 



= \n2 dt 



t^ 



1 + 



1 D 
8 ~ 24 



1 — imat — -m^a^t^ + 



(4.5) 



= 41n2 {-2m) ( ^ - ^m^a^ ) = Am ln2 (m^a^ - 1/12) 



2 2cosh(V2) 
^ ' t [2sinh(t/2)]^ 



n=l 



/ 



dt 



t^ 



1 + e 



1 £) 
8 ~ 24 



+ 



1 — imat — -rr?o?t^ + 



(4.6) 



1 i / 5 
imat\n2 - -m'^a'^t'^{2\n2 - 1) - -m^a^t^ 41n2 - - 

2 ^ ^ 6 V 2 



= 4 (-27r?) l/2m2a2 = -Am m'^a'^ 

And now calculate I2 on the contour that consists of a straight part from —ioo to — ie on the branch of 
the In function continuous with the left side of the cut, then a circle around the origin with radius e, and 
then another straight part from —ie to —ioo. First calculate the contribution of the straight parts. Since 
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the direction of integration on the two sides of the cut is opposite, only the contribution of the discontinuity 
across the cut remains, which can be expressed as 



2s J X ^ ' [-2isin(a;/2)]^ 



1 - x' 



1 D 
8 ~ 24 



+ ... 



-1 
2^ 



= / ^^2; ^ {-2m) i 

e 

oo oo 

- [ dx ^ {-ma) e"™"^ - Sni f dx ^ 



(4.7) 
(4.8) 
(4.9) 



Swi — (—ma) e 
2x ^ ' 



OO 

+ I dx Ani - {m^c? - 1/12) e"™"^ (4.10) 
J X 



+ 



OO 

- OO /. 

4wi Iny {m^a^ - 1/12) e"^ + 47ri / dy Iny {m^a^ - 1/12) e" 



(4.11) 



Considering that J dylnye ^ = r'(l), we get 




= Ani (2 + ema) e-"*"' - 47rj In(ema) {rn^a^ - 1/12) e" 

ema 

+ Ani {m^o? - 1/12) ^ r'(l) - J dy Iny e""^ 



(4.12) 



Our procedure is to take e ^ before we let ma go to infinity. This means that we can expand this 
result in e and discard all terms with positive exponents. 

= 47ri (2 + ema) { 1 - ema + 1/2 ^rr^c? ) - Ani In(ema) {m^a^ - 1/12) 

+ Ani {m'^a^ - 1/12) T'{1) ^"^"^^^ 



= Ame~'^ - Swie-'^ma + Qmm^a^ - Am {m^a^ - 1/12) In(ema) + 47ri (m^a^ - 1/12) r'(l) (4.14) 

Now we must calculate the contribution of the circle around the origin. The radius of the circle is e, 
and t can be expressed as ee'^. This means that Int = Ine + 20, is equivalent to itdO, and the integration 
around the circle from —ie clockwise means that 9 runs from 3/27r to — 7r/2. So for this part of the integral 
we have 

1 dO I 2s UU^0\ '"".t^tt (4-15) 
2s J t \ 2 J [2sinh(i/2)]'D ^ 

= 2^ In ^ ~lde e-^-* - 2 [e d9 ^^^'.''t^fi e"^-* (4.16) 

2 J [2sinh(t/2)]'D J [2 sinh(t/2)]^ ^ ' 
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(4.17) 



2 



4i In I (-27r) - ^m^a^^ + 4 J d9 (^r^ - imar^ + ^ ~ + - ) (^-1^) 



Here we drop the terms that contain t to a positive power, because on the circle t = ee*^ for t, and those 
terms will vanish as e — > 0. Thus we have 



o 1 e / 1 1 2 2 

-Sttz in m a 

2 V24 2 



+ A e d9 



imae e + 



1 1 
24 ~ 2' 



2 2 

m a 



(4.19) 



-Sm In^ f - -^m^a^ ) + 4 
2 V24 2 ' 



-2i 



1 



ima — e "^e + ( ^ - -^m^a^^ 4" 
-i V24 2 I 2 



(4.20) 



-4 / de 



-2i -i 



-Sttz In^ (^""■^™^"^) ~ 47ri e"^ + 8m mae~^ + 47r^ ( 7^ ~ ■^"^^"^ ) (4-21) 



= -47ri + 87ri e-^ma - 27r^ (m^a^ - 1/12) + 47rz (m^a^ - 1/12) In- (4.22) 

To obtain /(s), wc assemble the contribution of the pole at the origin and the cut along the negative 
imaginary axis to C(s), to first order in s. 

= 1 _ ,Us-n^ { ^1 + /3+ /4+ /5} (4.23) 
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[s + s2(21n(/ia)-2r'(l) + i7r) + 0(5^) 



s-^ 2m {m^a^ - 1/12) 



+ 4iri - Siri e~^ma + 6iri m^a^ - Ani {m^a^ - 1/12) In(ema) + Airi (m^a^ - 1/12) r'(l) 

- 47ri + 87Ti e-'^ma - 2tP- {m^o? - 1/12) + 47ri {rn^o? - 1/12) ln| (4.24) 



+ 47ri In 2 (m^a^ - 1/12) 



47r irr?c? 



= -^. \8 + s2(21n(/ia)-2r'(l) + i7r) + 0(5^)] 

— ZTTl 

{ 27ri (m^a^ - 1/12) + rn^o? - Airi {m^a^ - 1/12) \n{ma) (4.25) 
+ A-Ki {m^a^ - 1/12) r'(l) - 2t:'^ {rn^a^ - 1/12) } 



^ - {2iTi (rn^a^ - 1/12) 



-2-Ki 

+ s[{2 \n{na) - 2T'{1) + m) 2m {m^a^ - 1/12) 

(4.26) 

+ 2m rn^o? - Am {rn^o? - 1/12) In(ma) + Am {rn^o? - 1/12) r'(l) - 277^ {m^o? - 1/12) ] 

+ 0(5^) } 

= - (m^a^ _ 1/12) + s\-n?o? + 2(771^02 - 1/12) ln(m//x) ] + 0(5^) (4.27) 

From this expression, keeping only the leading behavior in (ma), and dividing by A-ko? (the volume of 
2-D de Sitter space) to obtain the zeta-function density, we find 

C(0) = (4.28) 



and 

A-K \ ' w? 

For the Minkowski zeta-function density in 2-D, we have from eq. (1.7) 



C'(0) = -^(l + ln^l . (4.29) 
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which gives precisely the values for ^(0) and C'(0) that we have obtained above. 



V. Conclusions 



In this paper we have, starting with the definition of the C-function for a minimally coupled massive 

scalar field in D-dimcnsional dc Sitter space, given an integral representation for it which is wcU-dcfincd near 
s = 0. We have then deformed the contour of integration so as to expose those terms which give the leading 
behavior as the de Sitter radius a tends to infinity. 



In this limit, dc Sitter space becomes Minkowski space, and wc verify, specializing to the case D = 2 for 
simplicity, that ^(0) and C'(0) ^^^'^ to the expected Minkowski values. This is a property that is not at all 
evident from the original integral representation, and, as pointed out in the introduction, is a property that 
might have been lost in the process of analytic continuation. It is both reassuring and a non-trivial check of 
our computations that the Minkowski limit is in fact recovered after analytic continuation. 



With these results in hand, one can consider the addition of a symmetry-breaking potential for the 
scalar field. One will then be able to study questions related to field quantization in de Sitter space in the 
presence of interactions that, in Minkowski space, would lead to symmetry breaking effects. In particular, 
one will be able to study the way in which particular time slicings may lead to symmetry restoration whereas 
others may allow the symmetry to remain broken [1]. Investigations along these lines are in progress. 
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